Shape coexistence in Lead isotopes in the interacting boson model with Gogny energy 

density functional 
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We investigate the emergence and evolution of shape coexistence in the neutron-deficient Lead 
isotopes within the interacting boson model (IBM) plus configuration mixing with microscopic 
input based on the Gogny energy density functional (EDF). The microscopic potential energy sur- 
face obtained from the constrained self-consistent Hartree-Fock-Bogoliubov method employing the 
Gogny-DIM EDF is mapped onto the coherent-state expectation value of the configuration-mixing 
IBM Hamiltonian. In this way, the parameters of the IBM Hamiltonian are fixed for each of the 
three relevant configurations (spherical, prolate and oblate) associated to the mean field minima. 
Subsequent diagonalization of the Hamiltonian provides the excitation energy of the low-lying states 
and transition strengths among them. The model predictions for the + level energies and evolving 
shape coexistence in the considered Lead chain are consistent both with experiment and with the 
indications of the Gogny-EDF energy surfaces. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60. Jz 
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I. INTRODUCTION 

The atomic nucleus is a physical system that exhibits 
a rich variety of intrinsic geometrical shapes: spherical, 
prolate and oblate. The coexistence and evolution of the 
different intrinsic shapes has been a major theme of inter- 
est in nuclear structure physics. It has been investigated 
extensively from both theoretical and experimental sides 
in the past decades [D-Q- 

In some specific regions of the Nuclide Chart, the en- 
ergies of the three intrinsic geometrical shapes bunch to- 
gether leading to the spectacular coexistence of three + 
states (including the ground state) in even-even nuclei. 
Neutron-deficient Lead isotopes present a nice example 
of the shape coexistence phenomena Q: In the 186 ' 188 Pb 
nuclei, the presence of three low-lying + states and other 
additional experimental data strongly suggest the coex- 
istence of spherical, prolate and oblate shapes. In the 
context of the nuclear shell model , the emergence of 
low-lying excited + states is traced back to the proton 
particle-hole excitation across the Z = 82 closed shell. 
The residual interaction between protons and neutrons 
is enhanced due to this cross-shell excitation, resulting 
in the lowering of the excited + states. In the vicin- 
ity of the N = 104 mid-shell, the effect is strengthened 
and has a stronger impact on excitation energies. For 
the 186 Pb nucleus, the three lowest + states are within 
a range of 700 keV and the two intruder + levels have 
the lowest excitation energy among the members of the 
Pb chain. The first excited + state is interpreted as a 
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proton two-quasiparticle (7^9/2) 2 intruder configuration, 
while the second excited + state could be interpreted as 
a proton four-quasiparticlc (7r/ig/ 2 ) 4 intruder configura- 
tion. These + states correspond to oblate and prolate 
equilibrium shapes. 

More quantitative results using large-scale shell model 
calculations can only be obtained in lighter nuclei. How- 
ever, for heavy nuclei including the Lead isotopes the di- 
mension of the shell model configuration space becomes 
exceedingly large and a truncation strategy preserving 
the essential ingredients of the low-energy spectrum is 
required. The Interacting Boson Model (IBM) [lj| has 
been successfully used for describing the low-lying states 
of medium-heavy and heavy nuclei, and presents a severe 
truncation of the full shell- model space (Tll - [l3 [. In this 
case, the building blocks are s and d bosons, which reflect 
the collective J" = + and 2 + pairs of valence nuclcons, 
respectively (llT - [l3j . 

Within the IBM, the description of intruder + states is 
based on the model by Duval and Barrett [l4], HH • They 
proposed to mix the normal (Op-Oh) configuration, com- 
prised of N bosons, with intruder configurations com- 
prised of N + 2n (n > 1,2,.. .) bosons, which takes 
into account the 2n-particle-2n-hole excitation across the 
closed shell. In the case of Pb isotopes with three low- 
lying + levels, the model consists of three different 
Hamiltonians corresponding to Op-Oh, 2p-2h and Ap-Ah 
configurations. The idea of configuration mixing in the 
IBM framework has been applied to spectroscopic anal- 
yses [H, fl~6l - fT9| . al gebr aic features [2(J [2l[, and geome- 
try and phases |22l424| associated with the shape coexis- 
tence observed in the Lead and Mercury region. In these 
studies, the parameters for the configuration-mixing IBM 
Hamiltonian have been extracted from a fit to the exper- 
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imental spectra and transition rates. 

The different configurations of the shell model are re- 
lated in the mean field language to the minima of the cor- 
responding mean-field deformation energy surface. The 
self-consistent mean-field method using microscopic en- 
ergy density functionals (EDFs) currently provides an ac- 
curate and universal description of nuclear ground-state 
properties and low-energy collective excitations, includ- 
ing mass, density distributions, surface deformation, gi- 
ant resonance, etc. The most popular EDFs can be of 
zero-range Skyrme [25| . finite-range Gogny [26[ as well 
as several parameterizations of the relativistic mean-field 
(RMF) Lagrangian [27|, [H|. The qualities and instabili- 
ties of the self-consistent description of shape coexistence, 
based on a series of Skyrme interactions, were examined 
in On the other hand, the so called NLSC RMF 

paramctrization has been tailored to describe the pro- 
nounced shape coexistence in Pb, Hg and Pt isotopes 
[30| . The Nilsson-Strutinsky method has also been used 
to study the neutron-deficient Pb and Hg isotopes [3l| . 

At the mean-field level, however, important symme- 
tries of the system are spontaneously broken. Therefore, 
to describe the spectroscopic properties of a given nu- 
cleus, one needs a systematic treatment of the dynamical 
effects associated with the restoration of the broken sym- 
metries and fluctuations in the collective coordinates. It 
is then necessary to project the mean-field solutions onto 
states with good symmetry quantum numbers and mix 
the different configurations. Configuration mixing calcu- 
lations, in the spirit of the generator coordinate method 
(GCM) have been performed for both Lead and Mer- 
cury nuclei, based on Skyrme [32|,[33j] and Gogny [34T - [3r| 
EDFs. 

A sound approximation to the full GCM configura- 
tion mixing calculation is represented by the solution of 
a five-dimensional collective Hamiltonian. Both vibra- 
tional and rotational mass parameters are obtained, from 
mean-field calculations, as functions of the quadrupolc 
collective variables. The collective potential is then taken 
as the total energy surface resulting from the mean-field 
approximation from which, the zero-point energies as- 
sociated with the rotational and vibrational motions are 
subtracted [37], IH| . This method can be also used for the 
description of shape coexistence phenomena based on an 
arbitrary EDFs, e.g., using the Gogny-DlS functional for 
Hg isotopes [m . 

More recently a comprehensive way of deriving the pa- 
rameters of the IBM Hamiltonian has been introduced 
[40j . By mapping the potential energy surface, obtained 
within the constrained self-consistent mean-field method 
with a given EDF, onto the expectation value of the cor- 
responding IBM Hamiltonian, the energy spectra and 
electromagnetic transition rates have been computed. 
This method has been successfully applied to various 
shape phenomena, including vibrational and 7-unstable 
[4l| as well as rotational deformed 42 1 nuclei, prolate- 
oblate shape transitions |43l| and to the study of the fin- 
gerprints of triaxiality [441 ] . 



In this paper we extend the method of Ref . [40( to take 
into account configuration mixing within the IBM. We 
will show, how the parameters of the configuration mix- 
ing IBM Hamiltonian can be determined without a fit 
to the experiment by using the microscopic input pro- 
vided by mean field energy surfaces in an appropriate 
way. Using this method, we are able to describe the emer- 
gence and evolution of shape coexistence. Our method 
is applied to the neutron-deficient Pb isotopes since the 
existence of three minima in some of them represents a 
quite stringent test of the model. Moreover, they are well 
studied both experimentally and theoretically offering us 
the possibility to benchmark our method with other pro- 
posals. Concerning the mean-field calculation, we use 
the Gogny-DIM 45| functional that was originally fitted 
to binding energies and radii. It has also shown good 
spectroscopi c prop erties as already exemplified in previ- 
ous studies [4a - l48| where it has been shown that DIM 
keeps essentially the same predictive power as the stan- 
dard Gogny-DlS EDF g|. 

The paper is organized as follows: In Sec. [Hi a brief 
review of the configuration mixing within the IBM and 
the geometrical interpretation is given. The mapping of 
the microscopic PES to the IBM one with configuration 
mixing is described and the way to extract the IBM pa- 
rameters is discussed in Scc. lIIII In Scc. lIVl the results of 
the diagonalization of the IBM Hamiltonian including en- 
ergy level systematics, the detailed level scheme and the 
B(E2) transition strength values for specific nuclei and 
the evolution of the spectroscopic quadrupole moment in 
the considered Pb chain are presented. Finally, Section 
fVlis devoted to the conclusions and work perspectives. 



II. DESCRIPTION OF THE MODEL 

We start with the self-consistent constrained Hartrec- 
Fock-Bogoliubov (HFB) approximation based on the 
Gogny-DIM EDFs. As constraints we use the mass 
quadrupole moments associated with the quadrupole de- 
formation variables j3 and 7 of the geometrical collec- 
tive model [50( . For a given set of collective coordinate 
variables q = (/?, 7), HFB calculations are performed to 
obtain the potential energy surface (PES) given by the 
HFB total energy denoted as E HFB (/3,j) (for details the 
reader is referred to [5l|). Note that, in some studies 
dealing with the five-dimensional collective Hamiltonian 
obtained from EDF calculations [U HI] , the PES is re- 
garded as the total energy obtained after subtraction of 
the rotational and vibrational zero-point energies to the 
HFB energy In our model the PES is simply the HFB 
energy and no zero-point energy correction are consid- 
ered. A typical example of such PESs is shown in Fig. Q] 
where the Gogny-DIM PESs in the (/3, 7) plane are given 
for the nuclei 182 " 192 Pb. The Gogny-DIM EDF calcu- 
lation produces a remarkable triple minima in 184 - 190 pb 
nuclei, where each local minimum is well isolated from 
each other. In the considered Lead nuclei, the spherical 
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FIG. 1: (Color online) Contour plots of the PESs as a function of the deformation parameters /? and 7 are given for the 
182 ~ 192 Pb isotopes. The plots denoted by DIM correspond to the microscopic mean field calculation with the Gogny-DIM 
EDF. The ones denoted by Mapped correspond to the mapped PESs used in fitting the IBM parameters. The color scale ranges 
from zero (the mean field ground state) to 3 MeV. 



minimum is always present while a prolate minimum de- 
velops from 182 Pb to 186 Pb, so does the oblate one. The 
prolate minimum becomes less significant from 188 Pb to 
190 Pb and finally disappears in 192 Pb. Within our model 
the Gogny-DIM EDF PES for an individual nucleus is 
mapped onto the corresponding IBM PES (to be dis- 
cussed below), as shown in Fig. [T] in the panels to the 
right of the Gogny ones. 

Let us turn to the IBM description of shape coexis- 
tence. In the present study, we consider the proton- 
neutron version of the IBM (usually referred as IBM- 
2) [ll], [l3[ since it takes into account proton excita- 
tions more explicitly than the original version of the IBM 
(IBM-1), which does not distinguish between proton and 
neutron degrees of freedom. The IBM-2 comprises the 
neutron (proton) s v (s T ) and d v (d^) bosons, reflecting 
the neutron (proton) = + and 2 + collective pairs of 
valence nucleons [Ill4l3| . The number of neutron (pro- 
ton) bosons, denoted as N v (N n ), equals the number of 
neutron (proton) pairs outside the inert core. 

To describe a system consisting of three different in- 
trinsic shapes, the Hilbert space is expressed as a direct 
sum of the orthogonal subspaccs for the normal (Op-0/i) 
and the two intruder (2p-2/i and Ap-Ah) configurations 



flil [l5j . The Hamiltonian of the system is written as 

H = P Q H Po+P 2 (H 2 + A 2 )P 2 

+P 4 (H i + A i )P i + H™ ix + H% x (1) 

where the Hi (i = 0,2,4) represent the Hamiltonians 
for the ip-ih configurations associated with the different 
intrinsic shapes, H^ iK (H^ ix ) are the interaction terms 
mixing the Op-Oh (2p-2h) and the 2p-2h (Ap-Ah) sub- 
spaces. The operators Pi are projectors onto the ip-ih 
configuration spaces and finally the A,; (i = 2,4) pa- 
rameters represent the energies needed to excite protons 
across the Z = 82 shell, which will be detailed later. 
We employ the Hamiltonian Hi written as 

Hi = eih d + KiQ% n '* ■ Q$" A , (2) 

where the first term fid = ndx + ndv represents the 
d-boson number operator while the second one is the 
quadrupolc-quadrupolc interaction between proton and 
neutron bosons. The quadrupole operator is defined as 
Q* pi = s\d p + d\s p + x P A d ldp} {2) (P = 7r ^)- In this 
and Xp,i are parameters. The Hamiltonian 
in Eq. ([2]) is taken in its simplified form in order to re- 
duce the number of parameters that are not directly de- 
termined from the PES. It keeps, however, the essential 
aspects of a more general IBM-2 Hamiltonian. 
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The mixing interaction terms H ] 
fined as 



2 i 
mix 



H, 



i-2i 



U)\ 2 % (44 + 8„8„) +i^2 2l {d\ ■ d\ + ■ dn)l?>) 



(i = 2,4) are de- for the expectation value of the Hamiltonian Hi 

Ei(p,~/) = ; m V KiN v N v 



1 



-2 i 



and 



-2 i 



stand for the mixing strengths. 



where 

In a shell model picture, the proton 2p-2h excitation 
across the closed shell Z = 82 creates one particle and 
one hole pairs in the Z = 82 — 126 and the Z = 50 — 82 
major shells, respectively. Since the IBM normally does 
not distinguish between particle and hole states, the 
2n,-particlc-2n-holc configuration comprises 2n additional 
proton bosons, and hence the model contains N v neutron 
bosons and N„ + 2n proton bosons. For the considered 
182-I92pk nycie^ the doubly magic systems 164 Pb and 
208 Pb are assumed to be the inert cores. As a conse- 
quence, the proton boson numbers are N v =0, 2 and 4 
for regular, 2p-2h and 4p-4h configurations, respectively, 
while N v varies between 8 and 11. 

A given IBM Hamiltonian can be related to the ge- 
ometrical model by the coherent-state framework [52j. 
The coherent state |$) represents the intrinsic wave func- 
tion of the boson system, and is written, up to a normal- 
ization factor, as 



i*>= n ( 



|o> 



(4) 



where the coefficients a pp are given by a p o = /3 p cos7 P , 
a p ±i = and a p ±2 = -^(3 p sm r y p . Here the parame- 
ters ft p and 7 P represent the axially-symmetric and the 
triaxial deformations for neutrons (p = v) and pro- 
tons (p = 7r), respectively. For simplicity we assume 
ftv = ft it — ft>B and 7„ = 7^ = 75. The ft parameter 
for the IBM is proportional to the one in the geometri- 
cal model. The proportionality coefficient is significantly 
larger than one due to the difference in the size of the 
model spaces [52| . On the other hand, the 7 variable can 
be the same for the IBM and the geometrical model. The 
PES for the IBM system of interest is given analyti cally 
as an energy expectation value of the coherent state [52| ■ 
The geometrical interpretation of the configuration 
mixing IBM was provided by Frank et al. (22J. The co- 
herent state in Eq. Q for a single configuration should 
be extended to be a direct sum of the coherent state for 
each configuration. The PES for the configuration mixing 
IBM is obtained as the lowest eigenvalue of the following 
3x3 matrix [22J 

Ea{P,i) ^02 (/?) 
E{p,i)=\ n 02 ((3) E 2 (ft }1 ) + A 2 Q M (ft) 

24 (/3) E 4 (ft, 1 )+A 4 



(5) 



where the Ei(ft, 7) (i = 0, 2, 4) in the diagonal part stands 



P 2 



B.i 



(1 + /3I,*) 2 



\l TjiXv.i + X*a)Pb,% COS 37 + jXv,iX-K,iPB,i\ ( 6 ) 



Here ft % B = with Cpj being the proportional- 

ity coefficient of the ft variable defined for the different 
mean-field minima associated with each configuration ip- 
ih, and N w j denotes the proton boson number in the 
ip-ih configuration. The non-diagonal entries f2j_2i(/3) 
(i = 2,4) represent the expectation values of the mixing 
interactions H^~ 2 \ given as 



^--2i(/3) 



y/N„4N„,j-l) 
1 + 



( fa 



{<»[- 

"B,i 

1 + Pba-%Pba 



L- 2 )(i+£L) 



^ l Pl,) 



■(7) 



Each of the microscopic PES, presented in Fig. [TJ is 
mapped onto the corresponding IBM PES, i.e., the low- 
est eigenvalue of the matrix in Eq. ([5]). Since the three 
local minima are well separated from each other, a set 
of parameters for each configuration are determined in- 
dependently from each others. First, the Op-Oh config- 
uration is assigned to the mean-field minimum with the 
smallest deformation. Then the 2p-2h configuration is 
assigned to the minimum with second larger quadrupolc 
deformation. Likewise the 4p-Ah configuration is associ- 
ated with the minimum with the third larger quadrupolc 
deformation. For each configuration, the parameters e,, 
Ki, Xu,i, X-kA and Cp ti in Ej(ft,y) of Eq. © are deter- 
mined, using the method of Ref. [41j . so that the topolo- 
gies, i.e., curvatures in both ft and 7 directions, around 
the corresponding minima are reproduced. For 186 Pb, for 
instance, the Hamiltonians for Op-Oh, 2p-2h and ip-ih 
configurations are assigned to spherical (ft = 0), oblate 
(ft « —0.2) and prolate (ft « +0.3) minima, respectively. 
Since the number of proton bosons is zero for all the 
considered Pb nuclei, the second term in Eq. ([2]) vanishes, 
and the parameters Ko, Xv,o an d X-nfi can be set to zero. 
Therefore, in the present study, the Op-Oh configuration 
always represents a pure U(5) limit of the IBM [1(| . 

The Aj parameters in Eq. ([5|) are constants depend- 
ing on the nucleus and they are fixed so that the en- 
ergy difference between the mean-field spherical and in- 
truder configurations is reproduced. These energy differ- 
ences between mean-field minima are denoted as 8E; = 



HFB 



(ft i ■ V • ) - E nFB (B° 



mn>7min) with 3 : 



n' Tmin) 



E 

being the coordinates that give the minimum for each of 
the ip-ih configuration in the HFB PES. These quantities 
should be in reasonable agreement with the observed + 
excitation energies. 

However, the values of the Aj derived from the esti- 
mation above should not be used in the spectroscopic 
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calculations with the Hamiltonian of Eq. ([T}, i.e., dif- 
ferent values of Aj should be used in Eqs. ((T|) and (|5|). 
From the original definition, the Aj (i = 2 or 4) repre- 
sents the offset energy added to the eigenenergies of the 
ip-ih Hamiltonian so that its ground-state + energy ex- 
ceeds that of the normal configuration by an amount that 
is roughly equal to the observed excited + ener gy a nd 
hence to SEi. More explicitly (cf. Appendix C of |15|). 



^(0 + ) + A. t = E o (0+) + SE i 



(8) 



where _E,(0 + ) represents the lowest (ground-state) + 
eigenvalue of the ip-ih Hamiltonian in Eq. @ . Note that 
the amount of energy gained by the mixing between nor- 
mal and intruder configurations is much smaller than the 
typical range of Aj values and is considered negligible in 
this rough estimate. In the considered Pb isotopes, since 
there is no deformation-driving term in the Hamiltonian 
Hq, the Eq(0 + ) energy is always equals to zero for the 
Op-Oh configuration. The lowest + eigen-energy com- 
prises the energy gained through the deformation at the 
mean-field level (equivalent to the depth of the minimum 
in the PES) and the extra correlation energy arising from 
quantum effects beyond the mean field. The A,; values 
determined solely by looking at the PES, do not take 
into account this quantum correlation energy, and hence 
is too small to describe correct spectroscopic tendencies 
consistent with the indications of the microscopic PESs. 

Let us consider, for example, the nucleus 186 Pb. The 
A 2 value derived from the PES, to be used in Eq. (|5|), 
is 4.014 MeV. Nevertheless, with this value, the intruder 
+ state becomes ground state after the mixing. This 
is apparently not consistent with empirical facts as well 
as with the indication of the microscopic PES. Since the 
2p-2h configuration gives -E^O" 1 ") = —3.676 MeV, to re- 
produce 8E2 = 1.208 MeV the A2 value to be used in 
Eq. (P should amount to A 2 = 1.208 - (-3.676) = 4.884 
MeV. The difference between the two A 2 values (= 0.870 
MeV), identified as the quantum correlation energy that 
the 2p-2h configuration gains through the diagonaliza- 
tion, seems so sizable as to change the conclusion. There- 
fore, for the spectroscopic calculations with the Hamilto- 
nian of Eq. ([TJ , we propose to use the formula in Eq. ((5J) 
to take into account the necessary quantum correlation 
effects. Also the Aj in Eq. (JTJ can be related to the 
ones in Eq. flSJ) by replacing Ei(0 + ) in Eq. ||8J by the 
deformation energy Ei(^ in , 7^j n ), and vice versa. The 
uncertainty in the parameters relevant to the configura- 
tion mixing has also been pointed out in Ref.[22J, where 
the PES of the configuration mixing IBM-1 Hamiltonian 
for Lead nuclei was analyzed. Although the parameters 
of the Hamiltonian give a good description of the spec- 
troscopy, only two (spherical and prolate) minima remain 
after configuration mixing in the 186 Pb nucleus Ij22j . This 
result seems to support our finding that the Aj values to 
be used in spectroscopic calculations may not at the same 
time give the IBM mapped PES similar in topology to 
the mean- field PES. 

To perform a fully consistent mapping of Aj in the 



present framework, the addition of some interaction term 
between like neutron bosons, such as of n v Q v ■ Q v type, 
to the Op-Oh Hamiltonian may solve the problem. The 
reason is that such term drives deformation and provides 
the energy which could compensate for the quantum cor- 
relation energy the intruder configuration gains. In fact, 
if one tries to put n v Q v ■ Q v with the realistic interac- 
tion strength k v = —0.013 MeV in the mapped Hamil- 
tonian Hq in Eq. (|2|) for 186 Pb, the Op-Oh configuration 
gives E o (0 + ) = —0.870 MeV, which is exactly the same 
as the correlation energy gained in the 2p-2h configu- 
ration. Nevertheless, since the microscopic Gogny-DIM 
PES suggests purely spherical minimum for the normal 
configuration, it is practically not possible to determine 
the strength parameter for such additional interaction 
term. Another possible solution which could work out is 
to map the angular-momentum projected PES onto the 
corresponding IBM PES. This could represent an inter- 
esting work for the future which is out of the scope of the 
present paper. 

The non-diagonal matrix elements, f2j_ 2 j(/3) in 
Eq. {TJ , concern the barrier between the mean-field min- 
ima but are only minor as compared to the diagonal parts 
in Eq. ((6|). Therefore, the parameter oj 1 ~ 2z can be intro- 
duced only perturbatively and is determined so that the 
barrier height for two different minima in the microscopic 
PES is reproduced. For the sake of simplicity, we assume 

i — 2i i—2i A — 2i 



III. MAPPED IBM POTENTIAL ENERGY 
SURFACES AND DERIVED PARAMETERS 

The mapped IBM and the microscopic Gogny-DIM 
HFB PESs are plotted in Fig. Q] for the nuclei i82-i92 pb 
In the case of 192 Pb the HFB approximation suggests two 
minima, and therefore only the Op-Oh and 2p-2h configu- 
rations are mixed in this nucleus. The location, relative 
energy differences as well as the energy barriers between 
the coexisting minima in the microscopic PESs are re- 
produced rather well in the mapped IBM PESs. Note 
that, due to the limited number of bosons, the mapped 
PESs are generally flat along the oblate axis. Although 
very shallow triaxial minima at 7 w 10° arc displayed in 
the HFB PESs of i88,i90 pb; in thc mappec i IB M PESs 

such minima arc approximated by axial ones. As a re- 
sult, some deviations of the barrier heights between the 
oblate and prolate minima occur for these nuclei. In or- 
der to describe the detailed ener gy s ystematics of quasi-7 
band, a boson three-body term |44fl is required which is, 
however, out of the scope of the present work. 

The IBM parameters, derived for the considered iso- 
topes 182 - 192 p D , are displayed in Fig. [21 Consistent with 
the evolution of the topology in the PESs shown in Fig.[TJ 
no rapid change with mass number is observed in these 
parameters. The comparison between the e parameters, 
for a given nucleus, in Fig. [Ha) reveals that eo is the 
largest, 64 is thc smallest while thc e 2 value is always 
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in between them. On the other hand, as a function of 
the mass number, eo looks parabolic with respect to the 
mid-shell nucleus 186 Pb, while €2 and £4 remain almost 
constant. Let us stress, that these boson number depen- 
dencies are consistent with the earlier phenomenological 
(see Ref. and references are therein) and microscopic 
fl2l [53| IBM-2 studies on collective structural evolution. 
The parameter K2 is, in general, larger than K4 as the 
model space of the latter contains a larger number of 
bosons. As functions of the mass number the \ param- 
eters, shown in Figs. [2jc) and[2jd), also display a weak 
dependence. Nevertheless, the sign of \u is always op- 
posite to the one of x-n- Their sum \v + X-* is positive 
(negative) for the oblate (prolate) 2p-2h (Ap-Ah) shapes. 
The Cp value (i.e., the scale factor for the j3 variable) 
does not change too much. This parameter is determined 
from the position of the axial minimum and the curva- 
ture along the j3 axis. Actually, the location of each axial 
minimum in the HFB PES in Fig. [1] remains almost the 
same in the corresponding IBM PES. 

FiguresJ^f) andHJg) show the energy offset Aj, defined 
in Eq. (|8]), and the mixing parameters for the Hamilto- 
nian fc^ , respectively. The magnitudes of both ui 02 
and ui are notably larger than those used in some fit- 
ted calculations within the configuration- mixing IBM-1 
model [HOI: ~ 10 kcV and w 24 w 20-30 kcV in the 
latter studies, while we have obtained ui 02 « 50 keV and 
lj 24 ps 200 keV. The present uj 24 value, which is partic- 
ularly larger than the one derived from phenomenology, 
implies that our microscopic EDF approximation sug- 
gests a complex topology of the mean-field PESs in the 
studied Lead isotopes in 7 direction. Therefore, it may 
require a mixing between the two intruder configuration 



spaces stronger than estimated from the pure fitting cal- 
culations. In particular, the mixing between the regular 
and 2p-2h configurations seems to be quite large in the 
case of 188 Pb. 

The offset energy Aj, depicted in Fig. HJg), roughly 
amounts to 4 and 8 MeV for the 2p-2h and Ap-Ah con- 
figurations, respectively. These values are approximately 
twice as large as the ones obtained in the IBM-1 phe- 
nomenology [l8[. One sees from Eq. ([5]), that a larger 
Aj energy is needed when the + eigencnergy of the 
intruder configuration is sufficiently large in magnitude 
compared to the + energy of the normal configura- 
tion. The quadrupole-quadrupole interaction for the in- 
truder configuration appears to be stronger in the present 
mapped IBM system than it is in the IBM-1 phenomenol- 
ogy. The intruder configuration gains a large amount of 
energy, giving rise to remarkable differences between our 
Aj values and the phenomenological [It], HH results. In 
fact, the derived K2 and K4 values are larger in magni- 
tude than those extracted from the fit. This may be due 
to the fact that the microscopic Gogny-DIM calculation 
(see, Fig. [T]) provides a pronounced minimum. 



IV. SPECTROSCOPIC RESULTS 

Having determined all the parameters required by the 
IBM Hamiltonian in Eq. ([1]) for each individual nucleus, 
the energy spectra and transition rates are calculated by 
diagonalizing the Hamiltonian within the enlarged model 
space consisting of the direct sum of the Op-Oh, the 2p-2h 
and the Ap-Ah subspaces. The results shown below are 
obtained without any fit to the experimental data, but 
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only from the Gogny-DIM HFB approximation and the 
mapping procedure described above. 

We have performed a diagonalization of the mapped 
IBM-2 Hamiltonian in the so called boson m-schcme ba- 
sis. The eigenfunction for each excited state gives rise to 
various spectral observables. In particular, the E2 transi- 
tion rates and the spectroscopic quadrupolc moments are 
important quantities by which one can gauge the emer- 
gence and the evolution of the coexistence and competi- 
tion between different shapes in the considered isotopes. 
For the E2 operator f < E2 \ we use the boson quadrupolc 
operator Q p ^, where the same parameter Xp,i as the one 
used in diagonalization is used, based on the idea of Cas- 
ten and Warner for the IBM-1 case (H^]. Within the con- 
figuration mi xing IBM framework this E2 operator can 
be written as [Ilflil 

f^ = J2e P ^T^ (9) 

P,i 

where e Pj i represents the proton and neutron boson effec- 
tive charges for each configuration. For simplicity, these 
charges are assumed to be the same (i.e., e„ ; ; = e Wt i = ej). 
For the effective charges, we have adopted the values 
given in Ref. [l| (i.e., e = 0.110, e 2 = 0.140 and 
e± = 0.170 eb). The effective charge should, in prin- 
ciple, be determined by taking into account core polar- 
ization effects. Such an effect could be renormalized in 
the effective charges used here, while a fully microscopic 
derivation of the boson effective charge still represents an 
interesting open problem. With all this in mind, the re- 
duced E2 transition B(E2; J — > J') between states with 
spins J and J', can be written as 

B(E2;J^ J') = -L-\ {J '\\f^\\J)\\ (10) 

where | J) and | J'} represent the wave functions of the 
initial and the final states with angular momenta J and 
J', respectively. 

The spectroscopic quadrupole moment Q (-0 f° r the 
state with spin J is given by 

Q {s) (J) = ^( J j o J j) (J\\t {E2) \\J), (ii) 

where use is made of the well known Wigner's 3-j symbol 

A. Level-energy systematics 

Figure |3] displays the theoretical (a) and the experi- 
mental [56[ (b) low-lying spectra as functions of the mass 
number. In the nuclei 184 - 188 pb, the relative location of 
the O2 and 2^~ experimental levels is nicely reproduced. 
Our calculations reproduce the correct location in en- 
ergy for these first excited + states, with the 0^~ level 
coming down as we approach the midshell nucleus 186 Pb 



and becoming the lowest-energy one at 186 Pb or 188 Pb. 
Both prolate and oblate minima become lowest in en- 
ergy for these nuclei (see, Fig. [1]) and therefore the resid- 
ual quadrupole-quadrupole correlation between neutron 
bosons and the intruder proton bosons become maximal, 
giving rise to these notably low-lying excited + states. 

The comparison between our results and the few avail- 
able data for the excitation energy of the 0;j" states re- 
veals that our values overestimate the experimental ones. 
This could be due to the fact that in the considered iso- 
topes, the third lowest-energy minimum in the mean-field 
PESs appears higher than expected from the experimen- 
tal point of view and also because of the level repulsion. 
Note, that the parabolic behaviour of the 0^" levels with 
respect to mid-shell is in good agreement with the rel- 
ative location of the three minima in the Gogny-PESs 
(see, Fig. [IJ: the three minima are closest to each other 
around 186 Pb while the second and third minima become 
less pronounced and only the spherical one remains as we 
approach the closed shells. 

The present calculations also reproduce the parabolic 
tendency of states with angular momenta J > 2. The col- 
lectivity of the intruder configurations becomes stronger 
and, as a result, the intruder states with J > 2 become 
most compressed around the midshell. Nevertheless, the 
change in all the calculated energy levels, including the 
excited + ones, takes place faster as compared with the 
experimental trend. Let us also stress that, similar to 
the situation observed for the 0^" levels, the calculated 
non-yrast 2^~, 4^~, 6 \ and spectra are more stretched 
than the experimental ones. 



B. Structure of eigenfunctions 

To interpret the dominant component in the calculated 
excited states and the structure of the wave functions, 
we show in Table HI the overlap probabilities of the basis 
states and the eigenfunctions corresponding to the three 
lowest-excited + states for all the considered Lead iso- 
topes. In all the isotopes the 0^ state corresponds to the 
spherical ground state with a nearly 100 % dominance 



TABLE I: Fraction of each configuration in the lowest three 
+ states of the considered 182-192 p D isotopes (in %). 



r 


Configurations 




184p b 


18bp b 


188p b 


iyu pb 


iya pb 




0p-0/i 


100.0 


99.8 


99.7 


98.6 


99.7 


99.6 


01 


2p-2h 


0.0 


0.2 


0.2 


1.3 


0.3 


0.4 




4p-4h 


0.0 


0.0 


0.1 


0.1 


0.0 






0p-0/i 


0.0 


0.1 


0.1 


1.4 


0.5 


0.5 


o+ 


2p-2h 


9.0 


14.6 


24.3 


65.5 


92.8 


99.5 




4p-4h 


91.0 


85.3 


75.6 


33.1 


6.7 






Op-Oh 


34.1 


0.9 


0.5 


0.7 


98.7 


93.2 


ot 


2p-2h 


41.2 


67.9 


67.5 


36.1 


1.2 


6.8 




4p-4h 


24.7 


31.2 


32.0 


63.2 


0.1 
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FIG. 3: (Color online) Level-energy systematics for 182_192 pb isotopes with mass number. Theoretical level energies resulting 
from the mapped IBM-2 Hamiltonian (a) are compared with the experimental (b) energies. The experimental data are taken 
from the ENSDF database [5f|. To guide the eye, each point has been connected. Solid, dashed and dotted lines stand for the 
lowest two excited + states, yrast states with J > 2 (2f, 4+, 6^ and 8*) and non-yrast states with J > 2 (2^, 4J, 62" and 
8%), respectively. Note that the experimental 2^ and 0^ excitation energies for 1S6 Pb ( 1S8 Pb) are 662 (724) keV and 655 (725) 
keV, respectively. 



of the Op-Oh configuration. For the nuclei I82,i84,i86pi^ 
the first excited + state is comprised predominantly of 
the Ap-Ah configuration, which corresponds to the pro- 
late minimum in Fig. Q] The extent of mixing between 
the 2p-2h and the Ap-Ah configurations for the first ex- 
cited + state becomes gradually stronger from 182 . 184 ph, 
to 186 Pb, which correlates well with the finding in Fig. Q] 
that the oblate minimum becomes more significant from 

182,184p b tQ 186p b _ 

Experimentally both 186 > 188 pb are regarded as the 
most spectacular examples of shape coexistence in the 
Pb isotopic chain. In this case, one sees a stronger mix- 
ing between different configurations in the first and the 
second excited + states. The 0^ state in 186 Pb is more 
or less clearly of Ap-Ah character while the two intruder 
configurations are mixed for th e Pit " state in 188 Pb. The 
earlier IBM-1 fitting calculation [17[ suggested almost the 
same predominance of the 0^" and the Ojj" eigenfunctions 
while the three configurations appear to be more strongly 
mixed for 186 Pb. On the other hand, the present re- 
sults for 186 Pb seem to be consistent with the ones ob- 
tained within the symmetry-projected GCM approxima- 
tion based on both the Skyrme-SLy6 [32| and Gogny-DlS 
[35| EDFs. In such studies [H, [35[, the collective wave 
function for the Oj" (O3") excited state is peaked on the 
prolate (oblate) side. We also find, that our results for 
the nucleus 188 Pb agree qualitatively well with the ones 
of previous symmetry-projected GCM studies [33|, Hf| 
where collective wave functions strongly peaked at the 
oblate and prolate sides have also been predicted. For 



the nuclei 190 < 192 Pb, there is almost no mixing between 
the different configurations for the three + states. In 
fact, the fraction of the Ap-Ah configuration is too small 
for them. 



C. Level scheme: 186 188 Pb nuclei 

In this section, we discuss in more detail the results 
obtained for the isotopes 186 Pb and 188 Pb which are the 
most distinct cases of shape coexistence in the consid- 
ered chain. We compare in Figs. @] and [5] our theoretical 
and the experimental energy levels and transition rates 
for these nuclei. The assignment of the calculated excited 
state to each band is done according to the predominance 
of a given configuration in the corresponding eigenstate 
and the E2 transition strength that exhibits a clear col- 
lectivity. 

For 186 Pb, in Fig. gj the calculated first excited 0^ 
state, predicted to be predominantly prolate, is quite 
close to the experimental value. From the experimen- 
tal point of view, such state has been identified Q as 
the oblate band head. On the other hand, our result in 
Fig. 0] is consistent with earlier predictions for the same 
nucleus within the symmetry-projected GCM approxi- 
mation based on the functionals Skyrme-SLy6 [32| and 
Gogny-DlS [35[. Actually, as seen from Table HI the Ap- 
Ah (prolate in the present IBM framework) component 
dominates 75.6 % of the 0^ state. A strong collective 
energy pattern is also predicted for this prolate band, 
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> 




FIG. 4: (Color online) Experimental and calculated energy 
spectra and B(E2) transition rates (in Weisskopf units) for 
the 186 Pb nucleus. Experimental energies and B(E2) values 
are taken from [SH, In the plot the experimental O3" and 
2f look nearly degenerated, but their excitation energies are 
650 keV and 662 keV, respectively. The experimental B(E2) 
of 510(120) and 6 (2) (W.u.) correspond to the 4f -> 2+ and 
2^ —¥ 0^ transitions, respectively. 



FIG. 5: (Color online) Same as Fig. g] but for the 188 Pb nu- 
cleus. Note that the theoretical 0^" and 2^ excitation energies 
are 1.086 MeV and 1.047 MeV, respectively. The following 
theoretical _B(E2) in the right panel are listed here to help 
identify the corresponding transition in the plot: -B(E2; 4^ — > 



2 J) = 236, B(E2; 2j -> 2 J) = 63, B(E2; 2^ 
B(E2;2+ -> 0+) = 1 (in W.u.). 



42 and 



with the ratio AE i+ /AE 2 + = 2.75. The B(E2) transi- 
tions among the members of this prolate band exhibit a 
collective behavior while the 2 + — > + E2 transition is 
very weak in the spherical band. Concerning the oblate 
band, the theoretical 0% excitation energy overestimates 
the experimental one. Note that, experimentally, this 0-^ 
state is recognized as the prolate bandhead 0. The ex- 
perimental 2^~ and 0% levels look nearly degenerated, and 
so does the present calculation except that the 2jj~ level 
lies slightly below the Ojj" level since the mixing between 
the two intruder configurations may be too strong. 

One notices from Fig. [5l that our model provides a 
similar level of quality in the description of the isotope 
188 Pb. Although the calculated excitation energy for the 
0^" state is a bit high, the calculated 2 state lies close 
to the experimental one. The present study also sug- 
gests, that the 0^" and the 0^" levels correspond to oblate 
and prolate configurations, respectively, which is consis- 
tent with symmetry-projected GCM calculations based 
on the Gogny-DlS EDF HH]. Nevertheless, the first and 
the second excited + states are experimentally in- 
terpreted as the prolate and the oblate bandheads, re- 
spectively. Moreover, the present study suggests a pro- 
nounced collective pattern for both the prolate (Ap-Ah) 
and the oblate (2p-2h) bands, and supports the experi- 
mental evidence for the strong E2 transition pattern in 
the band comprised of 2f , 4+, 6f, and 8f states. In our 
calculations, the two intruder + levels are rather close 
in energy, compared to the case of 186 Pb. In fact, among 
all the considered nuclei, the prolate-oblate energy dif- 
ference obtained from the Gogny-DIM PESs in Fig. Q] 



reaches its lowest value for 188 Pb. Due to the level re- 
pulsion, however, the excitation energy of the O3 state is 
larger than the energy difference between the spherical 
and prolate minima of the corresponding HFB PES in 

Fig.m 



D. Spectroscopic quadrupole moment 

The spectroscopic quadrupole moment Q^, computed 
according to Eq. (fTTj) . is shown in Fig. [5] func- 
tion of the mass number for the three lowest 2 + excited 
states of the considered Lead isotopes. In the case of 
182 Pb, Q (s) (2+) w cb and Q (s) (2+) « Q (s) (2+) w -2 
eb, reflecting the spherical and prolate character of the 
2^ and 2^,2^ states, respectively. The microscopic and 
the mapped PESs for this nucleus (sec, Fig. [1]) suggest 
a global spherical minimum and a well developed pro- 
late deformation. The two non-yrast 2 + states should 

originate from such a pronounced prolate minimum. For 
I84,i86p b; both the HFB and the mapped PESs Fig _[i] 

indicate the development of triple coexistence. The trend 
of the considered quadrupole moment changes accord- 
ingly. The Q (s) ( 2 f) (<2 (s) ( 2 J)) vamc is nearly -2 (+1.5) 
cb, suggesting that this state is prolate (oblate). From 
Fig. [5J one realizes that the prolate band consisting of 
the 2, 2^, 4*, 61 and 8^ states comes down in energy. 
On the other hand, our calculations suggest that the 
third band in 186 Pb, comprised of the 2f, Of, 4^, 6^, 
and 8^ states, originates from the 2p-2h oblate config- 
uration. Note, that the quadrupole moment for the 2jj~ 
state is positive. The same arguments apply to the nu- 
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182 184 186 188 190 192 



Mass Number 

FIG. 6: Calculated spectroscopic quadmpole moments Q' s ' 
for the lowest three excited 2 + states of the considered Pb 
nuclei as functions of mass number. Solid, dashed, and dot- 
dashed lines connect the calculated Q^ s ' values for 2+ , 2 J and 
2^~ states, respectively. 



clcus 184 Pb. 

A change in the spectroscopic quadrupole moments 
from 186 Pb to 188 Pb is also apparent from Fig. [6l In 
188 Pb, Q {s) {2f) = 0.19 cb while Q [s) {2f) = 0.36 eb. 
On the other hand, the Q^(2^) value becomes negative 
(= —1.10 eb). As can be observed from the level scheme 
displayed in Fig.[5j the 2\ state consists exclusively of the 
regular (spherical) configuration. The band consisting of 
the 0-^,2 2, 4^,6^ and 8^~ states emerges with predomi- 
nant 2p-2h oblate character while the one comprised of 
the 0^", 2^, 4j, d>2 anQl 82" states emerges with Ap-Ah pro- 
late character. It should be noted, however, that the 
spectroscopic quadrupole moment <3^ s ^(2^) for the 2 J 
state, assigned to the oblate band due to its stronger E2 
transition to the 0^ state, is quite close to the Q^(2^) 
value. This is mainly due to the fact that the mixing be- 
tween the different configurations is too strong for these 
two 2+ states: for the 2+ (2+) state, 58 (41), 32 (38) 
and 10 (21) % of its eigenfunction is composed of spheri- 
cal Op-Oh, oblate 2p-2h and prolate Ap-Ah configurations, 
respectively. The value Q^- S \2t t ) = —1.10 eb reflects a 
more clear prolate character, as the three configurations 
are less strongly mixed in this 2~£ state: 2, 38 and 60 % of 
the eigenfunction is composed of Op-Oh, 2p-2h and Ap-Ah 
configurations, respectively. 

For both i90,i92p bj we obtain that Q ( - s) (2+) is close to 
zero so that the 2f state is supposed to to be of spherical 
character. Our result seems to support the fact that the 
2^~ state is composed predominantly of the 2p-2h oblate 
configuration. This result agrees well with the corre- 
sponding PESs, shown in Fig. [1] for which the oblate 
minimum lies much lower, compared to 182_188 Pb, than 
the prolate one. Note also that (2+) < for !90,i92p b 
implying, that the 2^" state is prolate. 



V. SUMMARY 

To summarize, the emergence and evolution of the 
shape coexistence in the neutron-deficient Lead isotopes 
have been investigated within the configuration mixing 
IBM model with parameters extracted solely from a map- 
ping of the mean-field PESs obtained with the Gogny- 
D1M EDF. The diagonalization of the IBM Hamiltonian 
provides energy levels as well as transition rates between 
the excited states. It is important to emphasize that, 
although the IBM configuration mixing model contains 
many parameters, they can be determined unambigu- 
ously by relating the IBM PES for each configuration 
to the corresponding mean-field deformation minimum in 
the microscopic PES. No additional adjustment to exper- 
imental data is required. A potential difficulty and un- 
certainty of the fully consistent mapping concerning the 
offset energy A has been addressed and possible remedies 
for it have been discussed. 

The considered Lead nuclei present the most spectac- 
ular example of the coexistence of spherical, oblate and 
prolate equilibrium shapes. The relative location of the 
three associated + states were reproduced. In one of the 
most stringent tests, the 186 Pb nucleus, the present cal- 
culation suggested that the 0^" and the 0jj~ states are pre- 
dominantly of prolate (Ap-Ah) and oblate (2p-2h) nature, 
respectively. For the 188 Pb nucleus, another typical ex- 
ample with more available experimental data to compare 
with, the present work predicts the oblate bandhead as 
the first excited + state and the prolate band as the sec- 
ond excited + state. The calculated E2 transition pat- 
tern, albeit the quantitative deviation of the inter-band 
transitions from the experimental data, provides indica- 
tions of strong collectivity for the relevant prolate and 
oblate shapes. The experimental level-energy systemat- 
ics is well reproduced by our calculations. The study 
of the prolate-oblate dynamics has been complemented 
by looking at the spectroscopic quadrupole moment. Its 
value for different configurations and nuclei is consistent 
with the implications of other quantities and the sugges- 
tions of the mean-field microscopic calculations. 

Using the proposed methodology, many new re- 
search directions concerning complex shape dynamics 
are opened up. A possible application would be to an- 
alyze neighboring isotopic chains, Mercury, Polonium 
and Platinum isotopes. In particular, the study in 
the Platinum isotopes will help to disentangle if the 
single-configuration is the appropriate picture to describe 
those isotopes (see, e.g., Ref. [l9|, EH and references are 
therein). Other mass regions, including neutron-deficient 
krypton, selenium and germanium isotopes, and neutron- 
rich krypton, strontium and zirconium isotopes, which 
are also known as regions of shape coexistence [I| would 
be a potential target. 

The predictive power endowed to the model by the 
microscopic input makes possible the application of the 
present methodology to the study of exotic nuclei like the 
ones that will be experimentally accessible in the near 
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future. 
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